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Abstract 


We present JWST observations of the Crab Nebula, the iconic remnant of the historical SN 1054. The observations 
include NIRCam and MIRI imaging mosaics plus MIRI/MRS spectra that probe two select locations within the 
ejecta filaments. We derive a high-resolution map of dust emission and show that the grains are concentrated in the 
innermost, high-density filaments. These dense filaments coincide with multiple synchrotron bays around the 
periphery of the Crab's pulsar wind nebula (PWN). We measure synchrotron spectral index changes in small-scale 
features within the PWN’s torus region, including the well-known knot and wisp structures. The index variations 
are consistent with Doppler boosting of emission from particles with a broken power-law distribution, providing 
the first direct evidence that the curvature in the particle injection spectrum is tied to the acceleration mechanism at 
the termination shock. We detect multiple nickel and iron lines in the ejecta filaments and use photoionization 
models to derive nickel-to-iron abundance ratios that are a factor of 3-8 higher than the solar ratio. We also find 
that the previously reported order-of-magnitude higher Ni/Fe values from optical data are consistent with the lower 
values from JWST when we reanalyze the optical emission using updated atomic data and account for local 
extinction from dust. We discuss the implications of our results for understanding the nature of the explosion that 
produced the Crab Nebula and conclude that the observational properties are most consistent with a low-mass Fe 
core-collapse supernova, even though an electron-capture explosion cannot be ruled out. 


Unified Astronomy Thesaurus concepts: Supernova remnants (1667); Pulsar wind nebulae (2215); Pulsars (1306); 


Supernovae (1668); Core-collapse supernovae (304) 


1. Introduction 


The Crab Nebula is the first astronomical object identified to 
have resulted from a supernova (SN; Mayall 1939), and since its 
discovery almost 300 yr ago, it remains one of the most studied 
objects in the sky (e.g., Davidson & Fesen 1985; Hester 2008). It 
serves as a benchmark for understanding neutron stars and 
pulsars (Ostriker & Gunn 1969), and its well-established physical 
properties provide a critical anchor for calibrating models of 
neutrino-powered explosions (e.g., Sukhbold et al. 2016). 

We know the Crab Nebula resulted from a core-collapse 
explosion because it left behind a rapidly rotating pulsar. The 
pulsar converts its spin-down energy into a wind of relativistic 
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particles that, in turn, produce a synchrotron-emitting pulsar wind 
nebula (PWN). The ejecta in the Crab Nebula form a network of 
complex filaments consisting of ejected material that has been 
swept up and photoionized by the PWN. Studies of filament 
abundances and photoionization calculations support a relatively 
low-mass progenitor (e.g., MacAlpine & Satterfield 2008). 

The most recent estimate for the total mass of gas in the 
filaments is around 7 Mo, consistent with a ~9 Mo progenitor 
star and with most of the ejected mass accounted for in the 
filaments (Owen & Barlow 2015). The typical ejecta velocities 
of —.1200kms ! (Temim et al. 2009) account for less than 
10° erg of kinetic energy, and yet the explosion was an order 
of magnitude more luminous than a normal Type II SN (see 
Smith 2013). Chevalier (1977) argued that there must be a 
more extended SN remnant (SNR) beyond the visible Crab that 
would account for the missing energy, but despite intensive 
searching, no signatures of such a component have been found 
(Frail et al. 1995; Seward et al. 2006). While Sollerman et al. 
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(c) 


Figure 1. The NIRCam (panel (a)) and MIRI (panel (b)) imaging mosaic footprints are overlaid in blue on the Spitzer 24 um image of the Crab Nebula. Panel (c) 
shows the footprints of the MIRI MRS observations (small yellow squares), along with the fields of view of the two MIRI imaging observations taken simultaneously 


with the MRS (blue). 


Table 1 
JWST Imaging Observation 


Filter Wave. FWHM Exp. Time Extinction Correction Dominant Emission Components 
(um) (arcsec) (s) 
NIRCam F162M 1.626 0.055 870 1.184 + 0.068 [Fe n] 1.644, [Si I] 1.645, synchrotron 
F480M 4.834 0.164 870 1.030 + 0.010 Synchrotron 
MIRI F560W 5.589 0.207 555 1.027 + 0.010 [Fe n] 5.34, synchrotron 
F770W 7.528 0.269 1099 1.030 + 0.010 [Ar II] 6.99, synchrotron 
F1130W 11.298 0.375 422 1.049 + 0.018 [Ni rr] 10.68, [Ni m] 11.00, [He 1] 11.24, [Ni IV] 11.73, synchrotron 
F1500W 14.932 0.420 1055 1.032 + 0.012 [Ne III] 15.555, synchrotron 
F1800W 17.875 0.591 111 1.041 + 0.015 [S m] 18.713, dust, synchrotron 
F2100W 20.563 0.674 344 1.040 + 0.014 [S m] 18.713, dust, synchrotron 
F2550W 25.147 0.803 1099 1.035 + 0.013 [O Iv] 25.89, [Fe N] 25.99, dust, synchrotron 


Note. Values for filter wavelengths and FWHM of the PSF have been adopted from JWST User Documentation (https:/ /jwst-docs.stsci.edu/). Derivation of the 


extinction correction factors is outlined in Section 3.1. 


(2000) did detect a C IV A1550 absorption feature in the far-UV 
spectrum of the Crab's pulsar, the inferred mass and velocity 
are only 0.3 Mo and 2500 km s™', respectively, accounting for 
an additional energy of only 1.5 x 10° erg. 

Historically, such a weak explosion was theorized to be 
produced by an electron-capture SN (ECSN) resulting from the 
collapse of a super asymptotic giant branch (AGB) star with a 
degenerate O-Ne-Mg core (Miyaji et al. 1980; Nomoto et al. 
1984; Nomoto 1987; Janka et al. 2008). An ECSN was thus 
proposed to be the most likely culprit to produce the observed 
properties of the Crab (e.g., Nomoto et al. 1982; Kitaura et al. 
2006; Smith 2013; Yang & Chevalier 2015). More recent works, 
however, show that the Crab’s properties are consistent with a 
wider range of models, including iron core-collapse explosions 
(e.g., Gessner & Janka 2018; Hitomi Collaboration et al. 2018; 
Stockinger et al. 2020). Precisely what type of star and SN 
explosion produced the Crab Nebula remains unknown. 

Here we present JWST (Gardner et al. 2023) imaging and 
spectroscopic observations of the Crab Nebula that offer a 
uniquely detailed view of its PWN, ejecta filaments, and dust 
formation. With its unprecedented sensitivity and spatial 
resolution at infrared wavelengths (Rigby et al. 2023), the 
JWST observations allow us to probe the spatial variations in 
the Crab’s synchrotron emission that reveal insights about the 
acceleration of particles injected by the pulsar, isolate emission 
from dust to determine where the grains formed, and more 


accurately measure the ejecta abundances that constrain the 
Crab’s origin. 

The Letter is organized as follows. In Section 2, we describe 
the observations and data reduction, followed by a description 
of the imaging mosaics and the method for producing various 
emission maps in Section 3. We discuss the large-scale 
morphology of the ejecta filaments in Section 4, the dust 
distribution map in Section 5, large- and small-scale synchro- 
tron emission in Section 6, and the determination of the nickel- 
to-iron abundance ratios using photoionization models in 
Section 7. In Section 8, we discuss the implications of our 
results for the SN explosion and progenitor type, and we 
summarize our conclusions in Section 9. 


2. Observations and Data Reduction 


The Crab Nebula was observed with JWST in General 
Observer Cycle 1 under Program 1714 (PI: Temim). The 
majority of the SNR was imaged with the NIRCam (Rieke et al. 
2023) and MIRI (Wright et al. 2023) instruments, while 
Medium-Resolution Spectrometer (MRS) spectroscopy (Argyr- 
iou et al. 2023) was acquired in two chosen locations centered 
on the ejecta filaments. The footprints of the imaging and 
spectroscopic observations are shown in Figure 1, and the 
imaging filter parameters are summarized in Table 1, including 
the central wavelengths of the filters and the full width at half- 
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maximum (FWHM) of the point-spread function (PSF). All the 
data were downloaded from the Mikulski Archive for Space 
Telescopes at the Space Telescope Science Institute and can be 
accessed via DOI:10.17909 /6264-w578. 


2.1. Imaging with NIRCam and MIRI 


The NIRCam observations were carried out on 2022 October 
31 using the F162M and F480M filters. The Crab was mapped 
by a 3 x 1 mosaic, shown in panel (a) of Figure 1, using the 
default row and column overlap region of 10%. The size of the 
complete mosaic is 7/0 x 6/3. The chosen observing config- 
uration was the FULL array mode utilizing the FULL 3TIGHT 
primary dither pattern and 3 STANDARD subpixel dithers, for 
a total number of nine dithers per mosaic tile position. We used 
a BRIGHT readout pattern with five groups per integration for 
each dither position, for a total exposure time of 870 s per tile. 

The MIRI observations were performed on 2023 February 24 
using the F560W, F1130W, F1800W, and F2100W filters. The 
imaging mosaic consists of 4x3 tile positions, shown in 
Figure 1(b), and has an approximate size of 5/0 x 5/6. It was 
carried out with a four-point extended source dither pattern, a 
1096 row and column overlap region, and the FASTRI readout 
pattern. The exposure parameters for the F560W, F1130W, and 
F1800W filters all include a single integration with 50, 38, and 10 
groups, respectively. The F2100W observations were carried out 
with two groups of 15 integrations. The corresponding exposure 
time for each tile position is listed in Table 1. A dedicated 
background observation was taken at a blank-sky position off of 
the SNR, centered on coordinates a(J2000.0) = 05°34™38:3095, 
6(J2000.0) = 4-21?55'53" 53 and using the same exposure, dither, 
and filter parameters as the source observations. 

Additionally, during the MIRI MRS observations (see next 
section), simultaneous imaging observations were carried out in 
adjacent fields of view with the positions and orientation set by 
the telescope roll angle on the date of each observation. The 
footprint of these simultaneous imaging observations is shown 
in Figure 1(c). We chose a different imaging filter for each 
MRS grating change, leading to a total of three filters: F770W, 
F1500W, and F2550W. The exposure time for each of the 
filters is comparable to the total MRS exposure time and listed 
in Table 1. 

All imaging observations were processed using the JWST 
calibration pipeline version 1.12.0, the calibration reference 
data system (CRDS) version 11.17.2, and the CRDS context 
file jwst 1130.pmap. We used the World Coordinate 
System alignment tool JHAT (Rest et al. 2023) to align the 
NIRCam images to Gaia DR2. The skymatch step was 
turned off in the calwebb_image2 step of the pipeline since 
the extended emission from the Crab Nebula fills the entire 
field of view for some tiles. We used the MIRI dedicated sky 
observations to construct background images for each filter by 
coadding and sigma clipping the four individual dithers. This 
step is particularly important for removing the thermal 
telescope background that dominates at wavelengths longer 
than 15 um (Rigby et al. 2023). The background images were 
then subtracted from each of the level 2 images before they 
were combined into the final mosaics in the calwebb - 
image3 step. The final pixel scales of the NIRCam F162M, 
NIRCam F480M, and MIRI images are 0/031, 07063, and 
0/11 pixel !, respectively. All of the mosaics are shown in 
Figure 2, and the additional two-tile coverage obtained from 
the simultaneous imaging is shown in Figure 3. 
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2.2. Medium-resolution Spectroscopy 


The MIRI MRS spectroscopic observations were carried out 
on 2023 March 17 using all four channels (1—4) and the three 
grating settings (SHORT, MEDIUM, and LONG) to cover the 
entire wavelength range from 5 to 28m split into 12 
subbands, one for each channel/band combination. The 
MIRI MRS observations at the two filament positions shown in 
Figure 1(c) used the four-point extended source dither pattern, 
with one group of 100 integrations for each dither in the 
FASTRI readout mode, for a total exposure time of 1110 s per 
channel/band combination. The same configuration was used 
to observe a background region located to the south of the Crab 
Nebula, except the four-point POINT SOURCE dither pattern 
was used. 

We reduced the MRS data using version 1.12.0 of the JWST 
Calibration Pipeline, with versions 11.17.4 and "jwst 1154. 
pmap" of the CRDS and CRDS context, respectively. We 
processed all raw level 1B files through calwebb detec- 
torl to produce level 2A rate detector images. We create 
"master" detector background files for each MRS channel/ 
band combination by median combining the rate images before 
subtracting these from the corresponding detector images for 
each of the filaments. We then used calwebb spec2 to 
produce calibrated level 2B detector images, which included 
the nondefault detector-level residual fringe correction. We 
constructed cubes for each of the 12 MRS subbands using 
Spec3Pipeline, which implements the cube-building 
algorithm described in Law et al. (2023). As described in 
Argyriou et al. (2023), each of the subbands has differing 
spatial dimensions, with band 4C having a larger field of view 
than 1A. This is illustrated in Figure 4. 

We extracted spectra from a fixed rectangular aperture that is 
slightly smaller in size than the band 1A field of view from all 
the MRS subband cubes (see Figure 4). These spectra were 
processed with the additional spectrum-level residual fringe 
correction available in the JWST calibration pipeline to account 
for any residual fringing, in particular the high-frequency 
fringes in channels 3 and 4, thought to originate in the MRS 
dichroies (Argyriou et al. 2023), which are not completely 
removed by the detector-level correction in calwebb_spec2. 

Observing objects with extremely bright emission lines such 
as SNRs can produce several detector effects leading to 
spurious or misleading features in extracted spectra. These 
include the “pull-up/pull-down” electronic cross-talk effect 
(Dicken et al. 2022), light scattering in the detector across 
slices and spectral channels, and persistence (Argyriou et al. 
2023). After careful visual inspection of the extracted spectra, 
we found that all three affect the Crab Nebula spectra to some 
degree. Persistence and scattered-light effects produced false, 
unidentified broad and narrow emission lines, which we 
flagged and ignored in our analysis. The pull-up/pull-down 
affected the continuum in spectral channels with saturated 
emission lines. This manifests as the broad "dips" next to the 
brightest neon lines in the position 1 spectrum in Figure 5 that 
then also appear in other channels (e.g., the dip near 24 um). 
There is currently no correction for the pull-up/pull-down 
effect. 


3. Imaging Mosaics 


The composite image of the Crab Nebula created from all of 
the NIRCam and MIRI imaging mosaics is shown in Figure 6. 
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.F162M . F480M iln F560W 


F1130W F1800W F2100W 


Figure 2. The NIRCam and MIRI imaging mosaics of the Crab Nebula with the dominant emission components listed in Table 1. The size of each mosaic is 
4/9 x 5/5, and they are oriented with north up. 


F1500W 
Ne Ill 


Ar Il 
Ne Ill 
O IV + dust 


F2550W 
O IV + dust 


Figure 3. MIRI two-tile mosaics composed of images observed simultaneously with the MRS and the three-color image that combines all three filters. The position of 
the footprint relative to the entire Crab Nebula is shown in Figure 1(c). The size of each of the two MIRI mosaic tiles is 1/2 x 1/9. 
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1A: [Fe II] 1B: [Ni Il] 1C: [Ar II] 


2A: [Ne VI] 
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2B: [Ar III] 2C: [S IV] 


I"CENOD qm ae 


3A: [Ne II] 3B: [Ne V] 


Lem ue Lov 


3C: [Ne III] 


4A: [S III] 


4B: [Ne V] 4C: [O IV] 


LN A 


Figure 4. The extraction aperture for position 1 (white rectangle) plotted on images of emission lines in each of the 12 MRS subband cubes. The subband labels and 
species are shown in the top left-hand corner of each panel. The white lines in the bottom right of each panel represent 1" to highlight the increasing size of the field of 
view from subbands 1A-4C. The decreased spatial resolution from channels 1—4 due to the diffraction limit is also evident. North is up, east is left. 
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Figure 5. MIRI MRS spectra from the two observed positions centered on the ejecta filament in the equatorial region (position 1 in orange) and the southern filament 
(position 2 in blue). The two positions are shown in Figure 1. The detector artifacts that affect the shape of the spectrum are discussed in Section 2.2. 


Figure 7 enlarges particular regions of the composite image. 
For reference, at a distance of 2 kpc (Trimble 1968), one MIRI 
pixel corresponds to a spatial size of ~0.001 pc. These figures 
highlight the incredible detail present in both the synchrotron 
emission and the emission filaments while also demonstrating 
the large-scale structure: an inner “cage” of bright filaments is 
energized by the bright synchrotron from the pulsar. One gets 
the impression that the synchrotron emission then expands 
outward through gaps in the cage, in some cases extending 
beyond the outside edge of the ejecta filaments. Because of the 
combination of filters in this image, the filamentary emission 
shows a range of ionization levels and samples a range of 


filament densities. To provide insight into this, the dominant 
emission components sampled by each of the imaging filters 
are listed in Table 1, and the MIRI imaging filter transmission 
curves are overlaid on the MRS spectra in Figure 5. Below, we 
discuss how at least some of this complexity can be resolved by 
processing the images further to separate various components. 

For our study, we have chosen the minimum number of 
filters required to achieve the goal of mapping out the 
synchrotron emission, determining the distribution of dust 
grains, and mapping out the Fe emission. All of the images 
contain a significant contribution from synchrotron emission, 
but most also include some emission from the ejecta filaments. 
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Figure 6. STScI and NASA press release image of the Crab Nebula containing all MIRI and NIRCam images. The composite image shows the [S III] emission in red- 
orange, [Fe II] emission in blue, dust emission in yellow-white and green, and synchrotron emission from the PWN in smoky white. The colors assigned to different 
filters are blue (F162M), light blue (F480M), cyan (F560W), green (F1130W), orange (F1800W), and red (F2100W). Credit: NASA, ESA, CSA, STScI, T. Temim, 


J. DePasquale. 


The NIRCam F162M image contains synchrotron emission 
from the PWN as well as the Fe II 1.644 um and Sil 1.645 um 
emission from the ejecta filaments. The NIRCam F480M image 
is dominated by synchrotron emission with a minimal 
contribution from the ejecta filaments. This image is shown in 
Figure 8, and it provides a high-resolution and high dynamic 
range view of the PWN structure at IR wavelengths. 

All of the MIRI images show synchrotron emission from the 
PWN, but each also samples different emission lines from the 
filaments. The dominant emission line in the F560W mosaic is 
the [Fe II] 5.34 um line. While the F1130W mosaic contains 
mostly synchrotron emission, there is a small contribution from 
the ejecta filaments traced primarily by three different nickel 
lines ([Ni II] 10.68 jum, [Ni m] 11.00 um, and [Ni IV] 11.73 jum) 
and the He111.24 um line. The ejecta filaments in the F1800W 
and F2100W mosaics include emission from dust and the 
[S m] 18.71 um line. 

The imaging regions observed simultaneously with the 
MRS, shown in Figure 3, are dominated by [Ar II] 6.99 wm and 
[Ne m] 15.56 um lines for the F770W and F1500W filters, 
respectively, and by [O IV] 25.89 um and [Fe 1] 25.99 um for 
the F2550W filter. Clearly, the higher-ionization [Ne III] and 
[OIV] line emission extends significantly further to the 


southwest compared with the [ArH] emission and is related 
to the [O m] "skin" noted in previous Hubble Space Telescope 
(HST) studies (Sankrit & Hester ) Comparing the 
morphology of emission in these three images highlights the 
effect of ionization and density variations, with the low- 
ionization [ArII] appearing clumpy and structured while the 
[OIV] emission appears fluffy and diffuse. Although the 
[Ne m] and [O IV] panels look similar in extent in the black- 
and-white panels, the color image shows that their relative 
intensity changes with position as the green filaments are 
bounded by red toward the lower right. 


3.1. Method for Producing Emission Maps 


Here we outline the procedure used to create maps that contain 
isolated emission from dust and individual emission lines. The 
background emission for the MIRI filters was subtracted using 
the off-source blank-sky observations, as outlined in Section 
For the NIRCam images, we measured the background surface 
brightness in a circular aperture with a 10" radius centered on 
a(12000.0) = 05^34"40:815, 6(12000.0) = 4-22?02/48"347 and 
subtracted the measured value from the mosaics. An extinction 
correction was applied to all the mosaics using Ay = 1.08 + 0.38 
(De Looze et al. J) and the G23 average extinction curve 


THE ASTROPHYSICAL JOURNAL LETTERS, 968:L18 (25pp), 2024 June 20 Temim et al. 


Figure 7. Enlarged regions on the composite image shown in Figure 6. The [S III] emission is shown in red-orange, [Fe II] emission in blue, dust emission in yellow- 
white and green, and synchrotron emission from the PWN in smoky white. 
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Figure 8. The NIRCam F480M image of the Crab Nebula shown here is dominated by synchrotron emission from the PWN and has only a minimal contribution from 
the ejecta filaments. The image is shown on a log scale to maximize the dynamic range. The size of the mosaicked image is 4/9 x 5/5, and the orientation is north up. 


from Gordon et al. (2023) using Ry=3.1. To calculate the 
extinction correction factors for each wavelength, we used the 
Python package dust_extinction (Gordon 2023). The 
resulting extinction correction factor for each filter is listed in 
Table 1. Before producing the various emission maps, the images 
were convolved to the spatial resolution of the lowest-resolution 
image in the set using convolution kernels produced from PSF 
models generated by WebbPSF version 1.2.1 and the Python 
package pypher (Boucaud et al. 2016). 


3.1.1. Synchrotron Maps 


Our first step was to produce maps of the synchrotron 
spectral index and its uncertainty from the synchrotron- 
dominated F480M and F1130W images and use them to 
produce maps of the predicted synchrotron emission for each of 
the other MIRI filters. This process is outlined in Section 3.1 of 
Temim et al. (2012). We have assumed the currently quoted 
NIRCam photometric calibration uncertainties of ~1% for the 
F162M filter, 4% for the F480M filter,'^ and 5% for each of the 
MIRI filters (Dicken et al. 2024). The derived synchrotron map 
was then subtracted from the F560W, F1800W, and F2100W 
mosaics, leaving mostly emission from the filaments. We note 


8 https: //jwst-docs.stsci.edu /jwst-data-calibration-considerations /jwst- 
calibration-uncertainties 


here that since the NIRCam and MIRI data were acquired 4 
months apart, temporal changes in the synchrotron wisps 
(discussed in Section 6.3) produced artifacts in the FA80W— 
F1130W spectral index map in the central regions around the 
pulsar. These temporal changes did not affect regions farther 
from the pulsar and had a minimal effect on the emission from 
the ejecta filaments. We note that to examine the spectral index 
variations in the regions around the pulsar, we also produce a 
spectral index map using the contemporaneous MIRI F560W 
and F1130W images (see Section 6.2). 


3.1.2. Line and Dust Maps 


As summarized in Table 1, once the synchrotron component 
is subtracted from the remaining MIRI images, the residual 
emission in the F560W mosaic is dominated by the [Fe II] line 
emission and the F1800W and F2100W mosaics by dust and 
[S m] line emission. To produce separate emission maps for 
dust and [S M], we utilized the fact that the filter transmission 
curves for the F1800W and F2100W images conveniently 
overlap across a wavelength region that includes the 
[S m] 18.713 um line. This is demonstrated in the top panel of 
Figure 9, which shows the two filter transmission curves 
overlaid onto the MRS spectrum from position 2 and the 
relative contributions of various emission components shaded 
in different colors. The synchrotron component in blue has 
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Figure 9. The top panel shows the position 2 MRS spectrum with the F1800W 
and F2100W filter transmission curves overlaid and the different emission 
components shaded in different colors. The brightest lines contributing to the 
two filters are labeled. The bar chart in the bottom panel shows the contribution 
of the corresponding emission components to the total surface brightness in the 
F1800W and F2100W filters within the MRS field of view at position 2 (first 
two columns). It also shows the contribution of emission components to the 
weighted residual of the two filters (third column), where Residual — 
F2100W — 0.65 x F1800W. The dust map (fourth column) is given by 
subtracting the synchrotron component from the residual (Residual — 
Synchrotron), and the [S rif] image (fifth column) is given by [S I] = 
F2100W — Synchrotron — 1.65 x Dust. See Section 3.1 for details. 


been extrapolated from the F1130W band using the F480M-to- 
F1130W spectral index. The remaining continuum emission is 
assumed to arise from dust and is colored in red. The [S II] is 
shown in orange. The bar graph in the bottom panel of Figure 9 
depicts the contribution of each emission component to the 
total surface brightness within the MRS field of view at 
position 2. Even though the synchrotron component has been 
subtracted from the images, we show its contribution to the 
original image for reference. To isolate the dust emission, we 
subtracted the weighted F1800W image from the F2100W 
image, with the weighting factor of 0.65 chosen to completely 
remove the [S III] emission and isolate the dust emission. Other 
lines that contribute to the two images are colored green and 
labeled in the top panel of the figure. While their total 
contribution appears to partially cancel out in the bar chart, this 
would only be true if the spatial distribution of the emission 
lines was the same. Nevertheless, since the contribution of 
these lines to the total surface brightness is an order of 
magnitude lower than the contribution from dust, we are 
confident that the resulting residual image is dominated by dust 
emission. We note here that while the dust emission map 
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provides the spatial distribution of warm dust grains, the 
absolute values of the surface brightness are lower than they 
should be due to the dust emission in the F1800W being 
subtracted by the procedure we have used here. 

Our last step was to produce a map of the [S NI] emission by 
scaling the dust emission map (fourth bar in the bottom panel 
of Figure 9) up to match the total contribution of dust in the 
synchrotron-subtracted F2100W image and subtracting it from 
that image (rightmost bar in the bar chart). This new residual 
image is then dominated by [S II] emission with a small 
contribution from the other lines. Figure 10 shows the final 
resulting maps of the dust emission (middle panel) and the 
[S m] emission (right panel). The bottom panel of Figure 11 
shows the color-composite comparison of these two maps, 
where the dust emission is seen to be concentrated in the 
innermost filaments while the [S III] emission extends to larger 
radii. 


4. Morphology of the Filaments 


In this section, we present a brief assessment of the large- 
scale spatial distributions of the various maps constructed 
above in Section 3.1. Figure 11 shows color figures high- 
lighting various combinations of these maps. In all three panels, 
the [S III] image is used as a reference and is shown in green. 
The “cage” of inner filaments stands out from the more 
extended structure in all three of the comparisons shown. 
Interpretation of these comparisons is somewhat complicated 
because of the differing ionization stages of the various ions 
and the known correlation of ionization with density variations, 
especially in the bright filaments. Here we will concentrate on 
the larger-scale trends and leave more detailed smaller-scale 
filament comparisons to a future paper. 

The top panel in Figure 11 shows [O m] from HST (Loll 
et al. 2013) in red compared with [S m], so regions bright in 
both ions appear various shades of yellow. Because of the time 
separation between the existing HST data (circa 2000) and the 
JWST images, the proper motion of many filaments is obvious, 
and filament-by-filament comparisons are compromised. On 
larger scales, both ions display a similar filamentary morph- 
ology, but their spatial distributions are different. We see [O m] 
extending to greater radii, especially in the northeast and 
southwest regions aligned with the short axis of the nebula, 
while both ions extend along the major axis to comparable 
distances. The ionization potential for O^ to O** is 35.1 eV, 
while for S* to S‘*, it is 23.3eV. Hence, it is tempting to 
attribute the different spatial extents to asymmetrical ionization 
effects. However, the ionization balance calculations described 
in the next section reveal that [O mM] and [SII] emit over 
similar density and temperature ranges, so we cannot rule out 
the importance of varying elemental composition in contribut- 
ing to what is seen in this panel. 

The middle panel of Figure 11 shows the distribution of 
[Fe II] in blue to [S m] in green, both from JWST. With a low 
ionization potential of only 7.9 eV for Fe", this ion is expected 
to trace the densest cores of filaments that are being ionized 
from the outside inward by the synchrotron emission. Indeed, it 
is primarily the cage of bright ejecta that is prominent in blue in 
this panel, although it is interesting that fainter outer filaments 
along the major axis show up as well. The bottom panel of 
Figure 11 shows the warm dust map relative to [S m]. The dust 
emission primarily traces the bright inner filaments that likely 
have the highest densities. Along the major axis, some of the 
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Figure 10. The left panel shows the original F2100W image that contains synchrotron, [S m], and dust emission. The middle and right panels show the maps of warm 
dust and [S M] emission that have been derived from the synchrotron-subtracted F1800W and F2100W images (see Section 3.1). 


fainter dust emission seems to trace [Fe II] in the southeast, but Figure 12 shows a color-composite comparison of the JWST 
this behavior does not appear to be mirrored in the extended dust emission and the Herschel 70 um emission (PSF 
northwest region. The dust emission map is discussed in more FWHM = 5/6) that is likely dominated by the much cooler 
detail in Section 5 and the last paragraph of Section 6.1. grains. The brightest emission detected in the Herschel image 
also shows up in the new JWST data, highlighting that the 

5. Dust Emission warm dust is mainly distributed along the ejecta filaments that 


The frtevidence tor dustin the: Ciab Nebula came dronian host most of the cool dust material in the Crab. Explaining such 
observed IR excess above the spatially integrated synchrotron a wide dust temperature range may regue a scenarion which 
spectrum (e.g., Trimble 1977; Glaccum et al. 1982) and later the cooler dust is located deeper within the filaments, shielded 
from observed knots of optical extinction of the background from strong radiation, whereas the warmer dust may be located 
synchrotron radiation (Woltjer & Veron-Cetty 1987; Fesen & on the outskirts of the filaments, where they are heated. to 
Blair 1990). Estimates of the total dust mass in the Crab Nebula higher temperatures by the PWN's radiation. F igure 12 also 
based on the total integrated spectrum have varied consider- shows that the outermost filaments along the long axis of the 


ably, with a range of 0.02-0.5 Moe (Gomez et al. 2012; Temim Crab Nebula appear bluer and contain relatively more warm 
& Dwek 2013; Owen & Barlow 201 5). The total dust mass dust, while the pink, cooler dust emission peaks in the inner 
estimate was more recently refined to be 0.03—0.05 Mo after filaments. This could also be explained by more shielding of 
updates to the corrections for the Crab's synchrotron spectrum grains in the innermost massive and dense filaments and less 
and foreground interstellar dust emission (De Looze et al. shielding in the outer filaments with lower densities. — 
2019). The shielding of dust within dense filaments is consistent 
: : : 36 n 
Determining the spatial distribution of the emitting dust has with recent modeling efforts to explain the observed ""ArH 
been challenging; broadband imaging filters usually also emission in the Crab Nebula (Barlow et al. 2013) that require 
contain bright emission lines, and IR spectroscopy only exists high total hydrogen densities of 210 cm ` (Priestley et al. 
for select regions in the Crab. HST images showed that dust 2017) and 10 "cm ^ (Das et al. 2020) at the interface of 
extinction features are concentrated in the cores of the ejecta ionized and neutral regions within the Crab. Figure 12 also 
filaments (Blair et al. 1997; Sankrit et al. 1998), and spatially shows that the emission from the cooler grains peaks at the tips 
resolved Spitzer spectroscopy of select positions later con- of the Rayleigh-Taylor fingers, where the densities are higher 
firmed that dust emission does indeed reside in the filaments and more shielding from radiation is expected to occur. 
(Temim et al. 2012). Gomez et al. (2012) used Herschel Furthermore, the high densities in these regions may have 
imaging in broadband filters that is likely dominated by dust to stimulated efficient coagulation of grains, creating a size 
spatially associate a cool (28-34 K) dust component with the distribution that preferentially populates the large grain sizes 
ionized ejecta filaments that emit brightly in the optical (Fesen (a 2 0.1 um), consistent with polarized dust emission in the 
et al. 1992) and that host copious amounts of molecular Crab (Chastenet et al. 2022) and physical dust heating models 
hydrogen (Loh et al. 2012). (Temim & Dwek 2013; Owen & Barlow 2015; Priestley et al. 
The new JWST observations have allowed us to isolate the 2019). 
dust emission across the entire Crab Nebula, providing a dust The composition of dust in the Crab Nebula is still not well 
distribution map with unprecedented spatial resolution. This known. The three brightest ejecta filaments were detected in 
map is shown in the middle panel of Figure 10 and the bottom polarized emission with the SOFIA HAWC-+ instrument 
panel of Figure 11. The temperature of the dust grains emitting (Chastenet et al. 2022), and rather low polarization fractions 
in the MIRI F2100W image is approximately 55 K (Temim & (4%-10%) seem to indicate that at least part of the ejecta dust is 
Dwek 2013), so this map traces warm dust in an intricate composed of carbonaceous material. This is in line with 
network of ejecta filaments within the PWN. The method for observations of carbon-rich ejecta material found in the Crab 
producing this map is outlined in Section 3.1. (Gomez et al. 2012; Owen & Barlow 2015). One of the goals of 
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Figure 11. This figure shows color composites produced from the MIRI line 
emission images (see Section 3.1). The [O m] emission shown in the top panel 
was obtained in 1999-2000 with the WFPC2 instrument on board HST (Loll 
et al. 2013) and does not align on fine scales with the JWST images. See 
Section 4 for details. 


the MRS observations was to determine the composition of 
dust grains in the Crab Nebula. However, since the very bright 
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cooler dust 


Figure 12. A color-composite image with the warm dust component from 
JWST shown in blue and the Herschel 70 jum image dominated by cooler dust 
shown in pink. We note that the spatial resolution of the JWST/MIRI image of 
warm dust has a factor of 8 times higher spatial resolution than the Herschel/ 
PACS image (PSF FWHM of 0/674 vs. 5/6). 


emission lines in the Crab produced various still-uncharacter- 
ized detector effects that contaminate the continuum emission 
detected by MRS (see Section 2.2 and Figure 5), we leave a 
detailed study of the grain properties to a future paper. 


6. Synchrotron Emission 
6.1. Overall Morphology and the Synchrotron Bays 


The synchrotron-emission-dominated F480W NIRCam image 
of the Crab Nebula is shown in Figure 8. A remarkable level of 
detail in the PWN structure is evident in the image. The bright 
torus surrounds the pulsar in the center, with the outer nebula 
appearing to have a fibrous structure of intricate ripples and 
loops, seen in the enlarged images of select regions of the PWN 
in Figure 7, particularly panel 2. A long-known observational 
feature of the PWN is its elongation along the southeast/ 
northwest axis that is aligned with the pulsar’s jet. The long axis 
of the PWN measured ~7/6 in length, while the equatorial 
southwest/northeast axis that is aligned with the pulsar’s torus 
is ~5!5. 

This elongation of the Crab Nebula’s PWN was proposed to 
be caused by the pinching effect of the pulsar’s toroidal 
magnetic field (Begelman & Li 1992). However, some 
subsequent magnetohydrodynamic models for the evolution 
of the PWN showed that the elongated shape does not persist 
when the models are extended into three dimensions and that 
the close-to-uniform total pressure produces a more spherical 
PWN (Porth et al. 2014a, 2014b). In this case, the elongated 
shape would more likely be caused by an aspherical SN ejecta 
distribution or a disklike circumstellar material (CSM) 
distribution that confines the PWN in the equatorial plane, as 
proposed by Fesen et al. (1992). On the other hand, 3D models 
that assume an anisotropic distribution of the energy flux in the 
wind do produce an elongation in the PWN along the pulsar’s 
rotational axis (Olmi et al. 2016). 
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Figure 13. This image shows the F480M image of the Crab Nebula in blue that 
traces the PWN’s synchrotron emission and the dust emission map in red. The 
white arrows indicate the locations of indentations in the PWN that resemble 
the larger “bays” in the east and west equatorial regions that have been noted 
previously. 


Michel et al. (1991) noted an hourglass structure in polarized 
light images of the Crab Nebula and a highly organized “bay” 
structure in the east that was interpreted as scalloping in the 
outer nebula due to the PWN's interaction with the ejecta 
filaments. In this scenario, the filaments form a conducting cage 
around the PWN, and the wind expands outward between the 
dense filaments (Michel et al. 1991). Based on the measure- 
ments of the expansion of the synchrotron PWN, Bietenholz 
et al. (1991) and Bietenholz & Nugent (2015) also suggested 
that the relativistic gas is bursting through the net of filaments 
and accelerating into the surrounding low-density medium. 

These so-called synchrotron bays, the indentations in the 
synchrotron emission on the east and west sides of the PWN 
that give the nebula an hourglass shape, are shown in Figure 8 
and enlarged in panels 1 and 5 of Figure 7. The misalignment 
between the axis of the hourglass, which is oriented north— 
south, and the axis of the torus oriented in the southwest— 
northeast direction is a curiosity and has generated some debate 
on whether the global elongated morphology of the Crab 
Nebula is shaped by asymmetries intrinsic to the PWN, its 
confinement by asymmetric surrounding material, or some 
combination of both. 

The bays are proposed to be created through the interaction 
between the PWN and ejecta or CSM in the equatorial plane 
(Michel et al. 1991; Fesen et al. 1992; Li & Begelman 1992). 
Fesen et al. (1992) argued that the bays could be formed 
through the east-west confinement of the PWN by a disklike 
CSM structure, reminiscent of the ring in SN 1987A. The 
nebular magnetic field that wraps around this torus blocks the 
relativistic particles, causing the indentations in the PWN. The 
He-rich abundances of the east-west filaments that align with 
the bay structures (Uomoto & MacAlpine 1987; MacAlpine 
et al. 1989) offered additional support for this interpretation. 
Furthermore, a bipolar hourglass shape that is oriented north— 
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south is seen in the kinematics of emission lines from the 
filaments, such as [O m] A5007 emission (Smith 2003). The 
Crab's "chimney" feature also extends directly to the north of 
the inferred explosion center and is tilted from the southeast- 
northwest torus/jet axis (Davidson & Fesen 1985). 

Hester et al. (1995) pointed out that the east-west band of 
filaments may not be significantly influencing the structure of the 
PWN since other filament rings are also present in the nebula. 
They note that the east and west bays may look particularly 
prominent because the scalloping of the PWN by the filaments 
lies along our line of sight. There are, in fact, other baylike 
structures around the perimeter of the PWN. Figure 13 shows the 
F480M synchrotron emission in blue and the dust emission map 
(Section 5) overlaid in red. It can be seen that various other 
indentations or bays also coincide with the location of the densest 
filaments, as traced by the dust emission. The more prominent 
thicker filaments appear to be associated with larger indentations 
in the PWN and the thinner filaments in the southeast and 
northwest with smaller ones. This comparison suggests that the 
PWN is confined by multiple prominent filaments that lie in the 
wider northeast-southwest band in the plane of the torus and not 
only by the east-west running filaments, although the east 
equatorial filament does appear particularly prominent. The PWN 
is less confined in the northwest and southeast lobes that extend 
further out from the pulsar and show smaller indentations aligned 
with less prominent filaments. 

Recent 3D reconstruction of the structure of the Crab Nebula 
using imaging Fourier transform spectrometer data from the 
SITELLE instrument showed that the 3D volume follows a 
heart-shaped distribution that is symmetrical about the plane of 
the pulsar torus, with the two lobes separated by ~120° (Martin 
et al. 2021). This distribution is consistent with the confinement 
of the SN ejecta along the plane of the PWN torus. While this 
geometry may be produced by the expansion of a spherical 
PWN into CSM or ejecta that happen to be distributed in the 
plane of the torus, the alignment with the torus makes it more 
plausible that an axisymmetric nature of the PWN produces the 
asymmetry. 

Another example of a very elongated PWN in the SNR 
3C 58 (Reynolds & Aller 1988) shows the long axis to be 
aligned with the rotational axis of the pulsar (Slane et al. 2004), 
with optical ejecta filaments distributed in a roughly spherical 
distribution centered on the pulsar and extending to the radius 
of the PWN’s short axis (Fesen et al. 2008). 

It can be seen in Figure 13 that the synchrotron emission 
from the PWN extends well past the innermost dense filaments 
traced by the dust emission. 3D simulations of Blondin & 
Chevalier (2017) show that when the PWN expands into the 
outer steep part of the ejecta density profile, Rayleigh-Taylor 
instabilities can fragment the ejecta shell. The shocked pulsar 
wind breaks out of the shell and accelerates to larger radii in the 
freely expanding SN ejecta, leaving an inner shell of swept-up 
ejecta with about half its total mass enclosed within half of the 
PWN radius. The distribution of the dense filaments in 
Figure 13 appears consistent with this picture, particularly 
along the northwest direction, where the synchrotron emission 
extends well beyond the detected optical filaments (e.g., Loll 
et al. 2013). For the PWN to reach the steep part of the ejecta 
density profile, its deposited energy needs to exceed the SN 
explosion energy, consistent with a low-energy explosion for 
the SN that produced the Crab Nebula (Blondin & Cheva- 
lier 2017). However, this would require that about half of the 
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Figure 14. The top panel shows the synchrotron spectral index map computed 
from the coeval F560W and F1130W MIRI images. The map is shown on a 
linear scale with the index scale bar running along the bottom. The ejecta 
filaments have very flat indices and appear black in the image. The plot in the 
bottom panel shows the average surface brightness value within the white 
circular aperture shown in the top panel, a region chosen to avoid filament 
emission. The uncertainties on the values are dominated by photometric 
calibration uncertainties. The black line is a power-law model with a best-fit 
index of a = 0.47 + 0.10. 


kinetic energy we observe today in the Crab filaments is 
deposited by the PWN; thus, the original SN explosion was 
even lower, by about a factor of 2. Since a generous assessment 
of the current filament kinetic energy is <10°° (Smith 2013), 
this would require an extremely low SN explosion kinetic 
energy of only a few times 10° erg and would not account for 
the high observed luminosity of SN 1054. 

A similar filament morphology and PWN “blowout” will 
occur (Ostriker & Gunn 1971) if the PWN had encountered a 
thin shell produced by the interaction of the SN ejecta with 
significant CSM, as proposed by Smith (2013). During the 
CSM interaction phase, a large fraction of the SN kinetic 
energy is converted to radiation, making it even easier for the 
PWN energy to exceed the shell’s kinetic energy. This would 
also predict that the original SN explosion kinetic energy was 
actually somewhat more than the currently observed filament 
kinetic energy (i.e., more plausibly around ig^ erg), because 


13 


Temim et al. 


some of that explosion kinetic energy would have been lost to 
radiation to power the bright SN IIn event (Smith 2013). This 
scenario is also consistent with the lack of freely expanding SN 
ejecta and a blast wave outside the Crab filaments, since in this 
CSM interaction scenario, the forward shock has been 
decelerated by CSM. This latter scenario is discussed further in 
Section 8. 


6.2. Spectral Index Variations 


To explore the spatial variations in the synchrotron spectral 
index across the PWN, we computed a spectral index map 
between the coeval F560W and F1130W MIRI images (see 
Section 3.1). This map is shown in the top panel of Figure 14. 
Since the F560W image has a significant contribution from 
[Fe 1], the ejecta filaments stand out in the index map as black 
regions with low index values. The index map of the 
synchrotron emission itself is fairly smooth and seen to steepen 
with distance from the pulsar region, as previously seen in the 
optical (Veron-Cetty & Woltjer 1993) and IR (Temim et al. 
2006, 2012) observations. In this MIRI wavelength range, the 
index values are approximately 0.35 + 0.05 in the inner torus 
region to 0.80 + 0.07 in the outermost regions of the PWN. The 
quoted 1c uncertainties are dominated by the 596 photometric 
calibration uncertainty for the MIRI images and are 6 times 
lower than those measured for the spectral index map derived 
from the 3.6—4.5 um Spitzer images (Temim et al. 2006). 

The global spectral indices of the PWN measured at radio 
and optical wavelengths imply that a break in the spectrum 
occurs somewhere in the IR wavelength range (please see the 
broadband spectral energy distribution of the Crab Nebula in 
Figure 1 of Lyutikov et al. 2019). We searched for evidence of 
a break in the JWST data by selecting a region in the Crab 
Nebula that is dominated by the synchrotron component and 
relatively free of filament emission. This region was also 
chosen to be away from the inner torus region that exhibits 
temporal brightness changes on short timescales. It is shown as 
the white circle in the top panel of Figure 14. The average 
surface brightness values measured within this region for the 
NIRCam and MIRI images are plotted in the bottom panel of 
the figure. The black line is the best-fit power-law model with a 
3c uncertainty band shown in gray. All points are well fitted by 
a power-law model with a spectral index of 0.47 + 0.10, 
indicating that we cannot confirm a spectral break in the JWST 
wavelength range with the present calibration uncertainty 
values. 


6.3. Properties of the Wisps 


Highly structured continuum-emission features concentrated 
in the central regions of the Crab Nebula have long been 
understood to be associated with synchrotron radiation. Optical 
observations show that these features—ringlike “wisps” along 
with other more compact "knot" structures—show dynamical 
behavior on timescales of weeks to months (Scargle 1969), and 
simultaneous observations with HST and Chandra reveal 
corresponding structures and variations—both in position and 
brightness—in both the optical and X-ray bands (Hester et al. 
1996; Weisskopf et al. 2000), as well as in the radio 
(Bietenholz et al. 2001). These structures are understood as 
synchrotron radiation from the electron-positron wind as it 
enters the nebula from the termination shock, where the 
momentum flux of the wind is balanced by the pressure within 
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Figure 15. The first three panels are images of the inner pulsar region showing the structure of the synchrotron wisps, with a field-of-view size of 40" across and the 
scale in units of MJy sr !. The F560W and F1130W images were taken on the same day, while the F480M image was obtained 4 months before. The last panel shows 
the synchrotron spectral index map between the F560W and F1130W images. The spectral index values are labeled for several regions. While the absolute values of 
the spectral indices have an average 3c uncertainty on the order of ~0.15, the relative values in the index map have 3c uncertainties of only ~0.025. 


the nebula. The emission is highly polarized, and the magnetic 
field orientation follows the structure of the wisps. Theoretical 
models (e.g., Del Zanna et al. ) explain these wisps and 
knots as Doppler-boosted regions where the plasma bulk 
velocity points toward the observer. Wisps are typically located 
at intermediate latitudes, close to the equator, while the knot 
traces a polar outflow. 

Figure reveals these complex structures in the inner 
nebula as observed with JWST. The top left panel is from the 
NIRCam observation using the F480M filter and provides the 
highest-resolution image in the figure. Emission is seen from 
the pulsar itself along with a knot very close to the pulsar, 
labeled in the top right panel of Figure 15. The leading edges of 
the wisps in the northwest are Doppler-brightened due to the 
high-velocity equatorial flow and the geometry of the system, 
with the pulsar spin axis tilted into the plane of the sky by ~27° 
along a projected orientation roughly 126? north of east (Ng & 
Romani ). 
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The top right and bottom left panels in Figure 15 show MIRI 
observations of the central region obtained on the same day, 
taken through the F560W and F1130W filters, ~4 months after 
the NIRCam image. While the resolution of these images is 
lower, changes in the wisp structure are evident over the 
4 month period between these observations and the NIRCam 
observation. This is consistent with observations at other 
wavelengths that show wisp motions corresponding to 
velocities of —0.8—0.9c (Schweizer et al. ), typical values 
downstream of a shock in a magnetized wind with a ratio of 
magnetic to kinetic energy of o œ 1. 

The bottom right panel is a spectral index map produced 
from the F560W and F1130W maps (see Section ). The 
uncertainties on the relative spectral index values between 
different positions in the PWN are dominated by variations in 
the residual background emission and are significantly lower 
than those on the absolute values. The relative 1c uncertainty 
on the spectral index in individual pixels ranges from 0.0048 in 
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Figure 16. Left panel: spectral index of synchrotron emission due to pairs having an energy distribution given by a broken power law with lower energy index 1.35 
and high energy index 2.35 as a function of the break Lorentz factor ^, comoving magnetic field (orthogonal to the line of sight) Bj, and Doppler boosting term 
Dyoost (Where y is the bulk Lorentz factor of the emitting plasma and 5 is the component of the bulk velocity along the line of sight). The solid blue line is the 
analytical value at A = 8.3 jum; orange dots are values obtained from the ratio of integrated flux in the F560W and F1130W JWST filter bands. Right panel: difference 
in spectral index between a boosted (toward the observer) and unboosted emitting region as a function of the comoving magnetic field for two different values of the 
break Lorentz factor and two values of the bulk velocity of the boosted region at the effective wavelength of 8.3 um. 


the bright torus region to 0.024 in the outer PWN. This allows 
us to search for small variations in the spectral index across the 
Observed structures in the torus region. While the more diffuse 
regions of the nebula in Figure 15 have indices of ~0.35, there 
are significant variations for multiple structures. The region 
corresponding to the bright northwest wisp shows a consider- 
ably flatter spectrum, with o = 0.27, while that for the knot 
region is steeper, with a = 0.42. 

The cooling timescale for electrons emitting in the MIRI 
band (A~ 5 m) is fy; © 418. kyr. Thus, even if we 
assume a magnetic field of B ~ 400 uG, as expected down- 
stream of the termination shock in a highly magnetized (o œ 1) 
wind, the cooling timescale is much too long (given the typical 
timescales of the wisps) to explain the spectral index changes 
indicated in Figure 15. However, the broadband spectrum of 
the nebula indicates a break in the IR that is presumably 
associated with a break in the electron injection spectrum. The 
frequency of the observed synchrotron emission is 


(1) 


where Vsyn is the synchrotron frequency in the rest frame, ^y, is 
the Lorentz factor of the electrons, Boo is the magnetic field 
strength in units of 100 uG, and Dyoos = [TU — 8 cos0)]! 
is the Doppler boost factor associated with the pairs flow in the 
nebula (0 being the angle with respect to the line of sight). The 
observed frequency thus probes different parts of the electron 
spectrum depending upon the factor DboostB100. For example, 
for a uniform field strength, regions with a high Doppler factor 
probe lower-energy electrons than regions with unboosted 
emission. 

The flat spectrum of the northwest wisp, which is clearly 
Doppler-brightened, presumably probes lower electron ener- 
gies than for adjacent regions, thus sampling the electron 
spectrum below the break, while regions with steeper spectra 
probe electrons where the spectrum is steepening. The 
measured spectral indices of these features may be the first 
direct evidence that the spectral curvature in the injection 
spectrum of emitting pairs is tied to the acceleration process at 
the termination shock. 

Global spectral energy distribution modeling of the Crab 
Nebula (Bucciantini et al. 2011) indicates that emission comes 
from an underlying electron spectrum with a broken power-law 


2 
Vobs = DpoostUsyn x% 1.4Dy2Bioo Hz, 
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index distribution. The radio emission arises from particles 
with N (4,)) x 5; ?, while optical emission implies N (7) cx 
2 as (with ye being the Lorentz factor). Figure 16 shows the 
spectral index across the 5.6-11.3 jum wavelength range that is 
expected from synchrotron emission produced by a broken 
power-law particle distribution. The values of the spectral 
index o were derived using the analytical results of Gleeson 
et al. (1974). The orange points represent the index value 
calculated from the intensities in the F560W and F1130W 
filters, while the solid line represents the index values at an 
effective wavelength Aeff = 8.3 um, where yp, is the break 
Lorentz factor in the electron injection spectrum. 

Assuming that break energy m,c^^, is given by the ratio of 
the pulsar voltage and the pair multiplicity «, then 
Vee T. 9 x 10'°/k. Reasonable values of «25x 10°, 
^f. = 1.5 x 10°, and B œ 200 uG and observed wisp motions 
of ~0.9c (Schweizer et al. 2013) yield Aa=0.17 between 
boosted and unboosted regions, fully consistent with the 
variations seen in Figure 15. 

Emission from the knot region is much steeper than for the 
northwest wisp but is also Doppler-brightened. While the 
degeneracy between the effects of the Doppler factor, the 
magnetic field strength, and ^, makes it impossible to say with 
certainty, this is suggestive of a lower magnetic field from the 
emission regions of the knot and the large wisp to the southeast 
of the pulsar. This is consistent with pictures of the knot 
emission arising from high-latitude regions of the termination 
shock surface where the field is expected to be lower (Del 
Zanna et al. 2006). Another possibility is that different 
acceleration mechanisms are at play that produce different 
values of Y». For example, in the low-latitude regions where 
the wisps are thought to arise, the pulsar wind is expected to be 
striped and reconnection at the shock to play a major role in the 
acceleration. On the other hand, in the polar region associated 
with the knot, the pulsar wind should be unstriped and 
reconnection much less efficient (Lyubarsky 2003; Sironi & 
Spitkovsky 2011). 


7. Nickel-to-iron Abundance Ratios 


The existing estimated Ni/Fe ratios in select regions of 
the Crab Nebula range from ~6 times solar using near-IR 
lines to 75 times solar using optical lines (Henry 1984; 
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Figure 17. Left: the plot shows the [Fe II] 5.340 um (solid lines) and [Ni II] 6.636 um (dashed lines) line profiles as a function of the velocity shift with respect to the 
rest wavelength. The line profiles from position 1 and position 2 are shown in the green and orange colors, respectively. The solid gray curve shows the best-fit two- 
component Gaussian model to the [Ni II] 6.636 um at position 2, with the dotted and dashed gray curves showing the individual model components. The locations of 
the MRS pointings are shown in Figure |. Right: this figure shows line emission maps for the [Fe II] 5.340 uum and [Ni 11] 6.636 um lines extracted from the MIRI 
MRS data. The maps are produced by integrating the MRS cubes over the wavelength range that spans the detected line. While position 2 shows two distinct 
redshifted components (see left panel), position 1 contains distinct blue- and redshifted components with different spatial distributions. The white rectangles represent 
the spectral extraction regions used in the analysis. 


Table 2 
Tron and Nickel Emission Lines Used in the Analysis 
Position 1 (blue) Position 1 (red) Position 2 (red;) Position 2 (red) 
Line (jum) Cex Flux n;V Flux nV Flux n;V Flux n;V 
[Fe n] 5.3403 5.03e-09 61.3 : 1.5 4.48e48 14.6 + 1.1 1.07e48 334 2.41e48 214 x5 1.56e49 
[Fe n] 5.6739 2.19e-11 Ue ED EL EL 0.056 + 0.027 9.98e47 0.48 + 0.03 8.46e48 
[Fe 11] 6.7213 4.49e-10 3.76 + 0.14 3.87e48 0.96 + 0.10 9.88e47 1.16 + 0.25 1.19e48 14.2 + 0.2 1.47e49 
[Fe 1] 17.936 7.86e-09 30.4 + 1.2 4.7748 49+ 1.0 7.69e47 25+4 3.92e4 8 64 t4 1.00e49 
[Fe 1] 24.519 1.91e-09 7.2 X 0.5 6.36e48 0.66 + 0.13 5.83e47 6.0 + 0.9 5.30e48 16.2 + 0.6 1.43e49 
[Fe n] 25.988 3.95e-08 e EE Ue EL 11911 5.38e48 127 €7 5.75e48 
[Fe m] 22.9258 4.98e-08 25.3 + 0.7 1.06e48 5.9 t 0.5 2.47e47 8.5 + 1.0 3.55e47 28.3 + 1.1 1.18e48 
[Ni II] 6.636 2.42e-08 39.4 + 0.9 7.43e47 7.3 + 0.6 1.38e47 20.4 + 1.3 3.85e47 104 +2 1.97e48 


[Ni m] 11.002 5.87e-09 1.32 + 0.04 1.70e47 0.26 + 0.03 3.40e46 0.54 + 0.07 6.96e46 1.21 + 0.08 1.56e47 


Note. Fluxes at Earth are in units of 107! erg cm ? s^! and are converted to n;V by Flux x 2.1728e63/hv/C.x/n, to derive the total number of ions in the field 


of view. 


MacAlpine et al. 1989; Hudgins et al. 1990). Reconciling these [Fe 11] 5.340 um and [Ni 1] 6.636 um lines. It can be seen from 


early near-IR and optical estimates required either higher gas the plot that the line profiles of the Fe and Ni lines are the same, 
densities than have been generally accepted or the assumption confirming that these lines arise from the same velocity 
that the atomic data are inaccurate or incomplete (Hudgins et al. components. 

1990). The inconsistencies may also in part be due to the At position 1, two broad velocity components are evident 
complex ionization structure across the ejecta filaments (e.g., and are well fitted by individual Gaussian components, one 
Sankrit et al. 1998) that could not be properly accounted for in blueshifted by an average velocity of —500 kms! and the 
lower spatial resolution studies. The spectral and spatial other redshifted to ^ 750 kms '. The corresponding FWHM 
resolution and the sensitivity of the JWST MIRI MRS allow values of these two components are —260kms ! and 


us to measure these abundance ratios much more accurately. As 
will be seen below, the JWST data not only resolve the gas 
distribution that we are sampling but also provide many more 
Fe and Ni lines that make it possible to better constrain the 
densities and temperature of the emitting gas. 


7.140 km s, respectively. At position 2, almost all emission 
is redshifted, spanning a velocity range from 0 to 500kms_'. 
The line profiles at position 2 are well fitted by two broad, 
blended Gaussian components with centroid (FWHM) values 
of approximately 130kms ^! (200kms !) and 330km s`! 
(130 km s^ P). The spatial distributions for each of the two 


7.1. Spectral Line Measurements velocity components at both positions are shown in the right 
The MIRI MRS spectra seen in Figure 5 exhibit a wealth of panel of Figure 17, with the extraction aperture shown as the 
spectral lines arising from the filaments. In this initial analysis, white rectangle. It can be seen that the Fe and Ni lines have 
we restrict ourselves to a determination of the Ni/Fe abundance practically identical spatial distributions in addition to the same 
ratios, deferring more complete spectral line identification and velocity profiles, confirming that both lines arise from the same 
analysis to later work. As explained in Section 2.2, the MRS location. The other Fe and Ni lines exhibit the same properties, 
spectra in all four channels were extracted from a fixed aperture and their emission was fitted in the same way. The list of lines 
size shown in Figure 4. The left panel of Figure 17 shows the used in our photoionization models and their fitted flux values 
line flux profiles as a function of the velocity shift for the bright are listed in Table 2. 
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7.2. Photoionization Model Results 


We model the ionization balance following in part the 
approach in Laming & Temim (2020). We calculate the 
photoionization-recombination equilibrium for thermal gas 
residing in a Rayleigh-Taylor finger, constrained to have the 
same pressure as the surrounding PWN, taken to be 
1.9 x 10^? dyne cm ? (Porth et al. 2014b; Fraschetti & Pohl 
2017), which includes the contribution from electrons and 
positrons and the assumed magnetic field of 140 uG. 

The schematic geometry is shown in Figure 18. The PWN is 
taken to be a sphere of radius 1 pc, with a cone cut out of it 
where the thermal Rayleigh-Taylor finger gas resides. 
Referencing the lower left panel, we calculate the intensity of 
photoionizing radiation, J, in the Rayleigh-Taylor finger as 


ya & [^ re 2nd (cos Y) R?dR Q) 
An JR, Yoost; r?4- R? — 2rRcosq* 
where r?-— r? + R? — 2rRcosv, and in our case, 


cos 0; = —1 and cos 6, are given by the cone opening angle. 
The synchrotron emission per unit volume is e, and the PWN 
inner and outer radii are Rmin and Rmax respectively. 
Equation (1) becomes 


R max 2, 2. 
J= = inl ^ R 2Rr cos 0, RdR. (3) 
4r J Rmin r? + R? — 2Rrcos 62 
We put x = R — r cosÓ so that 
IG + R2 — 2Rr cos0)dR 
= f (x + r cos 0) In (x? + r? sin? 0)dx, (4) 


and with the help of Gradshteyn & Ryzhik (1965; Equations 
(2.733.1) and (2.733.2)), we arrive at 


J= " Uf Gas, 91) — f eas 92) — f (Rmin, 01) 


+ f (nin; 92)}, (5) 
where 
(R, 0) = 1R — r? cos20)In(R? + r? — 2rR cos) 
2 
= ler + r cos 0)? + 2r? cos? 6 
2 
+ 2r? sin 0 cos 0 arctan ( RIT = zl (6) 
sin 0 


It is useful to check that when 6, — 180° and 0, — 0°, 
Equations (A4) and (A5) of Laming & Temim (2020) are 
recovered. 

The PWN luminosity of 1.8 x 10° ergs ! in the energy 
range 300eV-10keV (Mori et al. 2004; Fraschetti & 
Pohl 2017) and radius of 1pc (assumed Rmin) give €= 
1.46 x 10? ergcm~*s~'. We calculate photoionization rates 
by using Equations (4) and (5) and integrating over a power- 
law spectrum with a photon index I’ = 2.1 (Mori et al. 2004). 
Collisional atomic rates are taken from Mazzotta et al. (1998), 
with updates to dielectronic recombination as given in Laming 
& Temim (2020). Photoionization cross sections are taken from 
the tabulation of Verner et al. (1996), with rates calculated by 
integrating over the PWN spectrum. Using an abundance set 
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for model V or VI in Owen & Barlow (2015), we iterate the 
calculation until the input electron density matches that at 
output. The sample ionization balance for Ni and Fe is given in 
Figure 18 and reflects a plasma at radius 0.9 pc in the PWN in a 
cone where cos 0 = 0.99, i.e., a filament radius of order 1% of 
the PWN radius, matching the geometry modeled in Sankrit 
et al. (1998). 

In Figure 18, the plasma particle density and temperature are 
coupled so that the gas pressure equals that of the PWN, given 
above as 1.9 x 10 ? dynecm ?. At high densities and low 
temperatures on the right-hand side of the right-hand panel of 
Figure 18, neutral atoms are found, with the equilibrium charge 
state increasing as the density decreases and temperature 
increases. The gray shaded region shows the approximate 
parameter space from which the emission we analyze arises. 
The high-density/low-temperature side is dictated by [Fe II] 
and [Ni II], to be discussed in more detail below, and the low- 
density/high-temperature side comes from the densities and 
temperatures associated with [O III] as observed by Fesen & 
Kirshner (1982). 

These limits are corroborated by an analysis of electron 
density diagnostic line ratios in [Fe II]. We calculate these using 
energy levels and A-coefficients from Deb & Hibbert (2011) 
and excitation rates from Bautista et al. (2015), including the 
lowest-lying 52 levels in our model. Figure 19 gives six line 
intensity ratios from [Fell] lines that are sensitive to the 
electron density, plotted for electron temperatures of 1350K 
(long dashed), 2400 K (thick solid), and 4270 K (short dashed), 
together with the observed ratios from positions 1 and 2, using 
the line fluxes in Table 2. At 1350 K and 4270 K, the 24.519/ 
25.988 and 17.936/25.988 ratios agree markedly less well with 
the other four, moving to higher/lower density as the 
temperature is decreased/increased. We use a temperature of 
2400 K and density of 3162 cm ^ to evaluate the [Fe n] and 
[Ni II] line intensities in Table 2 to give the abundance analysis 
in Table 3. For [Fem] and [Ni m], we take the same 
temperature and density, which is where these ions are 
maximized in Figure 18 (right panel). 

In addition to the [Fem] atomic data outlined above, for 
[Fe Mm] we use Zhang & Pradhan (1995) and Deb & Hibbert 
(2009), for [Nill] we use Cassidy et al. (2010, 2016), while 
[Ni m] is taken from Bautista (2001). Abundance ratios 
individually for Ni*/Fe* and Ni^'/Fe'" are given in 
Table 3, with each case corrected by the ionization balance 
to give a total element abundance ratio Ni/Fe. We find the Ni/ 
Fe number abundance ratio in the range of 0.156—0.277, 
corresponding to 2.8—5.2 times the solar ratio (now taken to be 
0.053; Scott et al. 2015), with the [Ni n]/[Fe I] generally 
registering higher Ni/Fe. The overlap between the density/ 
temperature ranges where the [Nill] and [Fell] lines are 
emitted is better than it is for [Ni III] and [Fe m]. [Ni III] comes 
from a slightly lower temperature/higher density than does 
[Fe m], and so a larger volume of plasma at the higher 
temperature would bias the Ni/Fe ratio downward. Never- 
theless, even the higher-end JWST-derived Ni/Fe abundance 
ratio of ~5 times solar is significantly lower than the previously 
reported ratios based on optical data. 

We compare our results with those from optical data of 
MacAlpine et al. (1989), who originally found the Ni/Fe 
abundance ratio to be of order 50 times the solar ratio. We have 
reprocessed their reported [Nin] 0.7378/[Fe m] 0.8617 flux 
ratios with our updated atomic physics models and find 
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Figure 18. Left: schematic diagrams (not to scale) of assumed PWN geometry. The thermal gas in a Rayleigh—Taylor finger resides in a cone cut out of the otherwise 
assumed spherical PWN. Synchrotron emission from the nebula irradiates thermal gas in the cone. The lower panel shows the geometry for the calculation. Right: 
ionization balance for Ni (top), Fe (middle), and variation of electron density and temperature (bottom) used in the Ni/Fe abundance ratio determination. 


Ni/Fe = 0.63, 0.50, and 0.57 for their positions a, c, and d. 
These amount to 9.4—11.9 times the solar ratio, the main 
difference being more modern and more complete atomic 
physics and our accounting for the ionization balance. Hudgins 
et al. (1990) measure the [Nill] 1.192/[Fe r1] 1.257 intensity 
ratio to infer the Ni/Fe abundance. Assuming an electron 
density «10? cm ?, they give Ni/Fe of about 0.3, or around 6 
times the solar ratio. However, the [Fe II] density diagnostic 
line ratio presented in their paper, the 1.294/1.257 ratio, 
indicates a density of 5 x 10? cm ?. Reevaluating this with the 
newer [Fell] model revises the density upward to around 
1 x 10* cm 7. At these densities, the inferred Ni /Fe abundance 
ratio becomes 30—60 times solar, because the [Ni II] 1.192 line 
is density sensitive. Hudgins et al. (1990) cast doubt on the 
reliability of this density diagnostic, and we would expect this 
higher-excitation line to be formed at a slightly higher 
temperature and lower density than the mid-IR lines. 

To have [Fell] appearing at a density of 10? cm ^, the 
ambient PWN pressure and/or the flux of ionizing radiation 
have to be decreased. The first does not seem plausible given 
the position within the nebula from which the Hudgins et al. 
(1990) spectrum was taken, but a reduction in ionizing flux 
could come about due to local extinction due to dust in the 
filaments. We find that a reduction by a factor of 0.1 moves the 
densities at which [Fe II] and [Nill] form to lower values by 
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about a factor of 0.5, so that the mid-IR lines observed by 
JWST are predicted to form at about 1500 cm ^? and the near- 
IR lines observed by Hudgins et al. (1990) at closer to 
1000 cm ?. At these lower densities, the JWST mid-IR lines 
give Ni-to-Fe abundance ratios in the range of 2.3-4.6 times 
the solar ratio (but with the higher end derived from [Ni II]/ 
[Fe N] being favored), the optical lines of MacAlpine et al. 
(1989) give 4.6-5.8, and the near-IR analysis of Hudgins et al. 
(1990) gives 8.2 times the solar ratio, bringing all of the 
estimates to a roughly consistent range. This flux reduction 
corresponds to 2.5 mag in the UV, which is consistent with V- 
band extinctions in the filaments reported by Grenman et al. 
(2017) and De Looze et al. (2019). Uncertainties on such 
abundance values have both observational and theoretical 
sources and are difficult to evaluate. We calculate the standard 
deviation of the number of Fe” ions in the field of view of each 
spectrum, taken from the four to six lines of [Fe II] observed, to 
find uncertainties of order 20%-50% for our JWST analysis. 
Similar uncertainties for the near-IR and optical results are 
plausible. 

We note here that if the Ni/Fe ratios vary spatially across the 
Crab's filaments, it is possible that the value we measure here 
may not be representative of the ejecta as a whole. However, 
based on the currently available data, there is no evidence that 
the Ni/Fe abundance ratios vary by more than ~50% across 
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Figure 19. Density diagnostic line ratios in [Fe II] at temperature 7 — 2400 K (solid), T — 1350 K (long dashed), and T — 4270 K (short dashed). Long and short 
dashed horizontal lines come from the red and blue velocity components of filament 1, and the dotted and solid horizontal lines come from the red1 and red2 velocity 
components of filament 2, respectively (see Table 2). All ratios indicate electron density in the range 1000-3000 cm ? for T = 2400 K. 


Table 3 
Derived Nickel-to-iron Ratios 


Derived Ratios 


Position 1 (blue) Position 1 (red) Position 2 (red;) Position 2 (red;) 
Line Ratio Cion bal. Ton Abund. x Solar Ton Abund. x Solar Ton Abund. x Solar Ton Abund. x Solar 
[Ni rr]/[Fe n] 1.61 0.152 0.245 4.6 0.162 0.261 4.9 0.147 0.237 4.5 0.172 0.277 52 
[Ni rr1]/ [Fe m] 1.13 0.160 0.18 3.4 0.138 0.156 2.9 0.196 0.221 4.2 0.132 0.149 2.8 


Note. The first column for each filament position gives the ion ratio, the second column gives the element mass abundance ratio that is corrected for the ionization 
balance, and the third column gives the Ni/Fe mass abundance ratio relative to the solar ratio. For reference, the solar value is Ni/Fe = 0.053. 


the Crab’s filaments. While the JWST data probe only two 8. The Nature of the SN Explosion 
select filament positions, we were able to isolate four distinct The Crab Nebula and SN 1054 have long been known to be 
ejecta velocity ranges that all have roughly the same Ni/Fe anomalous compared to what are usually considered “normal” 
abundance ratios. The other locations across the filaments that core-collapse SNe (CCSNe), and for decades (Nomoto et al. 
were observed in the near-IR and optical wavelengths have 1982) the Crab has been widely considered as the best observed 
Ni/Fe values that fall within a similar range (see above). The candidate for an ECSN explosion. This is based on a number of 
work of Sibley et al. (2016) indicates that the fractional Ni peculiar observed properties of the Crab: in addition to claims 
abundances relative to solar vary across the filaments (see their of an elevated Ni/Fe abundance ratio discussed above, these 
Figure 3), but it does not necessarily follow that the Ni/Fe include its very low explosion kinetic energy of less than 
ratios vary substantially. Additional spectroscopic observations 10°° erg, lack of any obvious a-element enrichment, moder- 
covering a range of filament positions will be required to ately low total mass in the filaments, low Fe abundance in the 
confirm a lack of significant variation in Ni/Fe across the filaments and low late-time luminosity from historical data 
filaments. (both of which imply a low ?Ni yield), location of almost 
In Section 8.1, we discuss the implications of the Ni/Fe ratio 200 pe out of the Galactic plane (implying an environment 
derived above. consistent with a relatively low-mass progenitor), and other 
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factors. However, a variety of theoretical work has continued to 
explore the end fates of stars that occupy the remarkable 
transitional range of initial mass (roughly 8-11 Mo) between 
those that die as planetary nebulae with white dwarfs and those 
that produce higher-energy CCSNe—ECSNe are not the only 
types of explosions that occupy this range. 

In the discussion below, we will refer to a few of these 
repeatedly, so it is convenient to define some acronyms that 
serve as shorthand for these explosive ends. We quote nominal 
single-star effective initial mass ranges for each, although these 
will obviously differ slightly from one model to the next 
depending on the adopted physics in stellar evolution and 
explosion models. 

ECSN. This is the ECSN discussed above, resulting from the 
collapse of a degenerate ONeMg core, with a stellar initial mass 
of roughly 8-9 Mẹ, and a low explosion kinetic energy of 
~10° erg. 

LMCCSN. This refers to low-mass Fe CCSNe that occur at 
the lowest end of the initial mass range immediately above 
ECSNe, at perhaps 9-10.5M,. These are also expected to 
produce explosions with low kinetic energy that in many ways 
resemble expectations for ECSN explosions. 

SDSN. This stands for Si deflagration SN, a special case of 
LMCCSN, discussed in detail by Woosley & Heger (2015). 
These progenitors have degenerate O and Si cores and ignite O 
and Si burning off-center in degenerate flashes as they 
approach Fe core collapse. Woosley & Heger (2015) find that 
for initial masses of 9.8—10.3 Mo, the Si deflagration flashes are 
robust and strong enough to eject much of the H envelope 
shortly before the Fe core collapse. Because of the pre-SN mass 
ejection, these events may show evidence of strong CSM 
interaction. 

CCSN. We will use this shorthand for Fe CCSNe that 
produce an explosion kinetic energy of 20.3 x 10?! erg (e.g., 
Sukhbold et al. 2016; Burrows et al. 2024a). These should 
occur above an initial mass of roughly 10.5-11 Mo, although 
the mass transition is not sudden and may differ from one set of 
models to the next. 

The first three above (ECSNe, LMCCSNe, SDSNe) are all 
expected to have a relatively low explosion kinetic energy of 
roughly 10? erg when the core collapses. Other explosions 
besides SDSNe may also produce strong pre-SN mass loss in 
degenerate flashes or super-AGB pulses or through wave- 
driven mass loss. Overall, there is much overlap in properties 
between some of these different models, which may make it 
very difficult to tell the difference observationally between a 
true ECSN and the low-energy explosions from lower-mass Fe 
core collapse. Below, we discuss a few of these factors that are 
particularly relevant to the Crab Nebula, including the Ni/Fe 
abundance ratio, explosion asymmetry and the neutron star 
kick, and possible CSM interaction. 


6.1. Ni/Fe Ratio 


A major observational goal of our study was to measure the 
Ni/Fe abundance ratios in the Crab Nebula filaments and to 
compare them to expectations from models of ECSNe and 
more typical Fe CCSNe. Previous studies, especially those 
using optical [Ni II] emission lines, had estimated high Ni/Fe 
abundance ratios in the Crab's filaments of 50—75 times the 
solar value (Hudgins et al. 1990). Theoretical predictions for 
nucleosynthetic yields of ECSNe have produced such high 
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Ni/Fe abundance ratios (e.g., Wanajo et al. 2009), seemingly 
favoring an ECSN model for the Crab. 

We investigated the Ni/Fe ratio in detail in Section 7, the 
main results of which were as follows. (1) The Ni/Fe 
abundance ratio in the Crab filaments estimated from our 
JWST spectra suggests a modestly elevated Ni/Fe ratio that is 
3-8 times the solar Ni/Fe ratio (assuming a 50% uncertainty). 
(2) We reanalyzed previous measurements including optical 
emission line strengths that led to very high estimates of the 
Ni/Fe ratio, and we found instead that with modern atomic data 
and reasonable ionization parameters, these observed optical 
line strengths are consistent with similarly modest Ni/Fe ratios. 
Using updated atomic data, the optical data give Ni/Fe ratios of 
9-12, but as shown in Section 7, accounting for dust in the 
filaments would lower these values. Overall, we suggest that 
observations of the Crab filaments provide a consistent estimate 
of the Ni/Fe abundance ratio that is on the order of 5 times the 
solar value, albeit with some mild differences in this ratio from 
one filament to the next. 

It is important to note that the expected or observed Ni/Fe 
ratio depends on the degree to which the products of explosive 
nucleosynthesis have been mixed with the ejecta and the 
unprocessed stellar envelope. One must be careful when 
comparing predicted nucleosynthetic yields for a specific layer 
or a specific burning region, a mixed emitting region, or the total 
yields averaged over all the ejecta. The filaments in the Crab most 
likely represent synthesized material that has been significantly 
mixed with the unprocessed outer envelope; the filaments contain 
several Mo of gas and show hydrogen emission (Sankrit et al. 
1998). Therefore, the Ni/Fe ratio produced by explosive 
nucleosynthesis in SN 1054 will be significantly diluted (and 
the resulting Ni/Fe ratio reduced) as it is mixed with the Ni and 
Fe in the unprocessed envelope that has a solar Ni/Fe ratio. Some 
theoretical model predictions for the Ni/Fe ratio refer to the pure 
results of the nucleosynthesis, and some refer to the total yield 
where these are mixed with the envelope. It is very important to 
distinguish these two when comparing observations to theory, 
and we refer to these as “unmixed” and “mixed” below. The 
relevant numbers for comparison with the Crab filaments will 
always be the “mixed” ratio. 

Very few studies report nucleosynthetic yields of stable Ni 
and Fe in ECSNe. Wanajo et al. (2009) examined the 
nucleosynthesis of ECSNe based on 1D models of an O-Mg- 
Ne core progenitor with an initial mass of 8.8 Mo (Nomoto et al. 
1984; Kitaura et al. 2006) and found Ni/Fe abundance ratios of 
24—57 times the solar value. These values are, however, 
unmixed. A more recent study of 3D Fe core-collapse explosions 
of zero and low-metallicity 8.1 and 9.6 Mẹ progenitors with 
tenuous mantles (that might at solar metallicity have O-Ne cores 
and explode as ECSNe) yielded Ni/Fe abundance ratios of ~20 
times solar (Wang & Burrows 2024a), though it is still unclear 
how similar the values would be for solar metallicity progenitors 
that explode as ECSNe. In all of the above models, the predicted 
ratio would be significantly lower if a stellar envelope is mixed 
in with the ejecta. For example, artificially assuming that a solar 
metallicity envelope is mixed with the ejecta of the "ECSN-like" 
models of Wang & Burrows (20242) reduces the Ni/Fe value to 
5-6 times solar, consistent with what we estimate for the Crab 
Nebula (A. Burrows 2024, private communication). While 
Wanajo et al. (2009) only provided the unmixed Ni/Fe values 
for their models, including the envelope would likely reduce 
their values in a similar way. Figure 20 shows the mixed and 
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Figure 20. The plot shows a comparison between the measured Ni/Fe ratios for the Crab Nebula and the predictions from some of the available explosion 
nucleosynthesis models. The green band represents the JWST-measured Ni/Fe ratios (relative to the solar ratio of 0.053 from Scott et al. 2015), including the 
estimated uncertainty of 50%. The near-IR measurements from Hudgins et al. (1990) and the optical ratios rederived in this work are shown as the gray bands. We note 
that the range of optical values does not include uncertainties in the measurements but represents a range of values measured at different positions in the Crab Nebula. 
These optical values are likely to be lower if dust in the filaments is accounted for (see Section 7 for details). The yields of the 3D models of Wang & Burrows (2024a) 
and Wang & Burrows (2024b) are shown as orange and blue solid vertical lines, respectively. The vertical lines reflect the variation in the expected ratio due to mixing 
with the stellar envelope (the dots are Ni/Fe values assuming no mixing with the envelope, and the crosses are values that assume the ejecta are mixed with the stellar 
envelope). The explosion energy (in units of Bethes) and resulting neutron star kick velocity are also labeled for each model (Burrows et al. 2024a; Wang & 
Burrows 2024a) and can be compared to the ~160 km s~' kick velocity of the Crab’s pulsar (Hester 2008; Kaplan et al. 2008). Nucleosynthetic yields from the 1D Fe 
core-collapse models of Sukhbold et al. (2016) are shown as blue diamonds for comparison. The orange dashed line connecting open orange circles represents the 
range of Ni/Fe ratios for the 1D ECSN models of Wanajo et al. (2009) that assume no mixing with the envelope. The orange dashed line connecting the orange pluses 
shows how these values might shift if we use a dilution factor of ~4 (similar to what is found for the ECSN-like models). For more details, please see Section 8.1. 


unmixed Ni/Fe predictions from these models, including the theoretical indications that this may happen more easily for 
estimated mixed values for the Wanajo et al. (2009) models that lower-mass stars (Woosley & Weaver 1995) and for asymmetric 
assume the same dilution factor as found for the Wang & explosions (Nagataki et al. 1997). While the more recent 1D 
Burrows (2024a) models. It can be seen from the figure that the explosion nucleosynthesis models of Sukhbold et al. (2016) do 
available ECSN yields in the mixed case are consistent with the not show this enhancement, a particularly asymmetric 3D model 
observed Ni/Fe values in the Crab. of an 11 Mẹ progenitor (Burrows et al. 2024a; Wang & 

The link between SN explosive nucleosynthesis and the Ni/Fe Burrows 2024b) does indicate a Ni/Fe enhancement of 
ratio in Fe CCSNe was investigated recently by Jerkstrand et al. ~4 times solar (envelope included), though the resulting neutron 


(2015b), following the inference of a high Ni/Fe ratio in the star kick velocity of —/700 km s^! is too high compared to the 
Type H-P SN 2012ec (Jerkstrand et al. 20152). The yield of stable 160kms ^! velocity of the Crab's pulsar (see Figure 20). 


nickel is dominated by the isotopes of ??Ni and 9?Ni, having The situation becomes more complicated if significant Fe 
neutron excesses of 0.034 and 0.033, respectively. The yield of and/or Ni contributions come from the innermost layers 
Fe is dominated by the synthesis of ?6Ni (which then decays to affected by neutrino weak interactions, thus changing the 
>Re through °°Co), with an initial neutron excess of zero. The neutron excess of the progenitor matter. These layers are too 
Ni/Fe ratio is therefore a diagnostic of the neutron richness of the small in mass (£102 Mo) to affect the total Ni and Fe yields 
layer that experiences explosive silicon burning. A solar ratio of for intermediate- and high-mass progenitors (i.e., CCSNe from 
Ni-to-Fe (0.053) corresponds to a neutron excess of 1.8 x 107°, 11 to 20M. progenitors) but may become relevant for low- 
close to the value found in the oxygen shell at solar metallicity. mass progenitors with lower yields of Ni and Fe that produce 
The next shell is the silicon shell, which has a neutron excess of a LMCCSNe. The recent work by Wang & Burrows (2024b) 
few times the solar value (Thielemann et al. 1990; Woosley & analyzing the nucleosynthesis of 3D Fe CCSN simulations 
Heger 2007; Jerkstrand et al. 2015b). The normal CCSNe that (Wang & Burrows 2023; Burrows et al. 2024a) found that the 
present Ni/Fe ratios of 3—5 may therefore be interpreted as stars 9 Mo progenitor model with some imposed initial velocity 
in which the dominant Fe and Ni emission lines arise in perturbations produced more neutron-rich material, leading to a 
condensations dominated by material that was explosively burned final Ni/Fe abundance ratio of ~8 times solar. This value 
and ejected from the deep-lying silicon-shell material. There are reflects the ratio for the entire ejecta, including the stellar 
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envelope, so it is appropriate to compare it with the observed 
ratio in the Crab’s filaments that might contain mixed-in outer 
ejecta layers or unprocessed CSM. The fact that two different 
9 Mo Fe core-collapse models of Wang & Burrows (2024b) 
result in very different Ni/Fe ratios (see Figure 20), with the 
model with imposed initial velocity perturbations producing 
higher Ni/Fe, suggests that a wide range of Ni/Fe ratios is 
possible for models in this mass range of progenitors and that 
these models are consistent with our reported Ni/Fe ratio for 
the Crab. 

The observationally derived Ni/Fe ratios for the Crab 
Nebula and the predictions from the models discussed above 
are summarized in Figure 20. It can be seen that the modestly 
elevated Ni/Fe ratio in the Crab is consistent with the “mixed” 
Ni/Fe ratios for either an LMCCSN or an ECSN. These SNe 
also provide low enough explosion energies (listed in units of 
Bethes (B) in Figure 20) to be consistent with the Crab 
(Burrows et al. 2024a; Wang & Burrows 2024a). The pulsar 
kick velocities resulting from the 3D models are also included 
in the plot. As discussed in Section 8.2, a stronger line of 
evidence for an Fe core-collapse explosion for the Crab may be 
the higher asymmetry of these explosions that can reproduce 
the measured kick velocity of the Crab’s pulsar (see Figure 20 
and Section 8.2). Other recent evidence includes the honey- 
comb-like structure revealed in a 3D reconstruction of the 
Crab’s ejecta (Martin et al. 2021) that may be produced by 
expanding plumes of radioactive ?9Ni-rich ejecta seen in Fe 
CCSN simulations (Stockinger et al. 2020). Some recent efforts 
focus on the analysis of emission lines from lighter elements 
that may have significant differences between ECSNe and 
LMCCSNe (Hosseinzadeh et al. 2018; Jerkstrand et al. 2018; 
Hiramatsu et al. 2021), which may be a useful avenue for 
further constraining the explosion that produced the Crab 
Nebula. 


6.2. Asymmetry and Pulsar Kick 


Simulations of ECSNe and Fe CCSN models predict 
different ejecta distributions in the resulting SNR. Fe core- 
collapse models predict higher asymmetries in the ejecta in the 
critical seconds after the explosion that can impart pulsar kick 
velocities exceeding the 160 kms ! velocity of the Crab's 
pulsar. The asymmetry in the postshock ejecta leads to an 
acceleration of the pulsar in the opposite direction from the 
heaviest ejecta species, including Fe (Wongwathanarat et al. 
2013; Janka 2017), an effect now observed in multiple SNRs 
(Holland-Ashford et al. 2017; Katsuda et al. 2018; Bhalerao 
et al. 2019). 

Recent 3D SN explosion simulations investigating neutron 
star kicks through the gravitational tugboat mechanism 
revealed a potential challenge to interpreting the Crab Nebula 
as an ECSN. The steep density gradient in the O-Ne-Mg core 
of ECSN progenitors leads to explosions that are much more 
symmetric than Fe CCSNe; thus, ECSNe are expected to 
produce low neutron star kick velocities of only a few to tens of 
kms ^! (Gessner & Janka 2018; Burrows et al. 2024b). This is 
much lower than the 160 km s^! measured for the Crab’s pulsar 
(Hester 2008; Kaplan et al. 2008). So it seems that the Crab's 
pulsar received a kick that is perhaps too fast for an ECSN but 
too slow for a higher-mass CCSN. Indeed, simulations of 
~9 Mo Fe CCSNe can achieve more asymmetric explosions 
with SN energies and kick velocities in the range of what is 
observed for the Crab Nebula (e.g., Stockinger et al. 2020; 
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Burrows et al. 2024a), as shown in Figure 20. The observed 
motion of the Crab's pulsar may therefore favor an LMCCSN, 
an Fe CCSN with a low progenitor mass. In order to still be 
consistent with an ECSN, the observed pulsar velocity would 
require an alternative kick mechanism, such as anisotropic 
neutrino emission or a breakup of a binary system. 


8.3. CSM Interaction 


If SN 1054 was not a true ECSN but was instead a low- 
energy Fe core-collapse event arising from the lowest initial 
masses to have this fate (LMCCSN), it is worth discussing how 
it might have produced such a high luminosity. These 
LMCCSNe are expected, like ECSNe, to produce relatively 
low-energy explosions (kinetic energy of ~10°° erg), they may 
be equally consistent with the observed Ni/Fe ratio, and they 
may be more consistent with the observed motion of the pulsar. 
With a progenitor resembling a super-AGB star or red 
supergiant, such explosions would presumably yield relatively 
low-luminosity SNe II-P. 

Indeed, both ECSNe and low-mass Fe core-collapse events 
have been invoked to explain several relatively low-luminosity 
SNe II-P, including events similar to SN 2005cs. This is, 
however, at odds with the fact that SN 1054 was extremely 
bright, being visible during the daytime for a few weeks with 
an estimated absolute visual magnitude of —18 mag (3 mag 
brighter than SN 2005cs and other faint SNe II-P, even though 
the Crab's kinetic energy is even lower than estimates for this 
class of extragalactic SN). Omand et al. (2024) explore the 
possibility that the luminosity of SN 1054 could be explained 
by a pulsar-driven SN model, similar to those used for 
superluminous SNe, but this scenario still leaves the question 
of the very low observed kinetic energy. Smith (2013) 
proposed that the low kinetic energy, high peak luminosity, 
thin filaments, and several other properties of SN 1054 and the 
Crab could be reconciled if the event was a low-energy 
explosion that was brightened significantly by strong and 
immediate interaction with CSM, likely resembling the 
subclass of events called SNe IIn-P (Mauerhan et al. 2013; 
Smith 2013). In this scenario, the strong CSM interaction 
occurring during the main SN event would cause much of the 
ejecta and CSM to be swept up into a very thin shell, which 
would cool radiatively into a very thin and dense layer and may 
even form dust in the postshock layers. This cold massive shell 
would later be fragmented by the expansion of the PWN, 
eventually producing the Crab filaments we see today. This 
model, however, invokes the ejection of a massive CSM shell 
shortly before the final Fe core collapse. 

In this context, the phenomena discussed by Woosley & 
Heger (2015) may be of interest. They found that stars with 
initial masses just above the ECSN range (corresponding to 
initial masses of 9-10.3 Mo in their models) suffer off-center 
degenerate Si-burning flashes. In some cases, these Si flashes 
may be powerful enough to cause the ejection of much of the H 
envelope shortly before Fe core collapse, leading to a bright 
CSM-interaction-powered SN. Above, we referred to these 
events as SDSNe, with "SD" referring to the Si deflagration 
that precedes Fe core collapse. Woosley & Heger (2015) 
estimate that this type of pre-SN mass ejection happens 
robustly in the initial mass range of 9.8-10.3 Mẹ in their 
models, although similar but less powerful flashes may occur 
over a wider mass range. Such a scenario may provide a viable 
alternative to an ECSN to explain the origin of the Crab Nebula 
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and makes specific predictions about the timing and energy of 
pre-SN mass-loss events. This scenario can be tested directly if 
light echoes from SN 1054 can be found because spectra of 
these echoes could confirm the prediction of a Type IIn-like 
spectrum (Smith 2013), probably similar to other SN IIn-P 
events. 


9. Conclusions 


We presented the analysis of JWST infrared observations of 
the Crab Nebula composed of MIRI and NIRCam imaging 
mosaics that sample iron, sulfur, dust, and the PWN's 
synchrotron emission, as well as smaller field-of-view imaging 
that samples argon, neon, and oxygen emission. The observa- 
tions also include two MIRI MRS pointings centered on the 
Crab’s inner Rayleigh-Taylor filaments, and they reveal spectra 
rich in emission lines, including many lines from iron and nickel. 
The unsurpassed sensitivity and spatial resolution of the JWST 
data reveal the Crab Nebula in unprecedented detail and have 
allowed us to study the large- and small-scale properties of the 
PWN's synchrotron emission, isolate emission from dust grains 
in the ejecta filaments, and determine nickel-to-iron abundance 
ratios in the ejecta to compare them with explosion and 
nucleosynthesis models of Fe CCSNe and ECSNe. Below, we 
summarize our main findings and conclusions. 


1. All emission lines from the ejecta filaments display the 
familiar filamentary morphology but with different spatial 
distributions that likely reflect a combination of ionization 
effects and varying elemental composition across the 
filaments. We used a differencing technique with the 
longer-wavelength MIRI images to produce a high- 
resolution map of dust emission across the Crab Nebula. 
We found that the emission from the mid-infrared-emitting 
warm grains is concentrated in the bright, highest-density 
inner filaments. A comparison of the dust emission map 
from JWST with a Herschel image sampling emission 
from cooler grains shows that the outermost filaments 
along the Crab's long axis have relatively warmer dust, 
possibly due to less shielding from radiation in the outer 
filaments. The emission from cooler grains dominates at 
the dense tips of the Rayleigh-Taylor fingers, reflecting 
either enhanced shielding or a grain size distribution that is 
weighted toward larger grains. 

2. The comparison of the synchrotron-dominated JWST 
images with the dust map that traces the densest filaments 
of ejecta also leads to interesting observations. Various 
indentations or "bays" are evident in the synchrotron- 
dominated JWST images and are seen to coincide with 
the location of the densest ejecta filaments, as traced by 
the dust emission map. The PWN appears to be confined 
by multiple prominent filaments distributed in a wide 
band oriented along the torus and less confined along the 
Crab's long axis, where thinner filaments and smaller 
PWN indentations are seen. The PWN extends well past 
the densest filaments along both axes, a morphology 
consistent with simulations of synchrotron “blowout” that 
occurs when the PWN energy exceeds the SN explosion 
energy (Blondin & Chevalier 2017) or when it encounters 
a thin shell produced by the interaction of the SN ejecta 
with significant CSM material (Ostriker & Gunn 1971; 
Smith 2013). 
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3. A synchrotron spectral index map produced from the 
F560W and F1130W MIRI images shows index values 
ranging from 0.35 + 0.05 in the inner torus to 0.80 + 0.07 
in the outer PWN. Flux densities extracted from a 
filament-free region of all NIRCam and MIRI images are 
well fitted by a power-law model with a spectral index of 
0.47 + 0.10. With the current MIRI calibration uncertain- 
ties, we were not able to confirm a break in the global 
synchrotron spectrum at JWST wavelengths. However, 
due to the very small uncertainties in the relative 
synchrotron spectral index across the PWN, we measured 
index changes across various small-scale structures in the 
PWN's torus region. The region corresponding to the 
bright northwest wisp shows a considerably flatter 
spectrum (a = 0.27) than the more diffuse regions of 
the torus (—0.35). This change in the index can be 
explained by Doppler boosting of emission from particles 
with a broken power-law distribution. In this scenario, the 
flat spectrum of the Doppler-brightened northwest wisp 
probes lower electron energies than for adjacent regions, 
thus sampling the electron spectrum below the break. The 
measured index change may be the first direct evidence 
that a break in the injected particle spectrum is tied to the 
acceleration process at the termination shock. The knot 
region and a large wisp southeast of the pulsar have a 
steeper spectrum (œ= 0.42) than the surroundings, 
suggesting either a lower magnetic field in these regions 
or a different acceleration mechanism. 

4. We detected several Ni and Fe lines in the MRS spectra 
of the ejecta filaments and used photoionization models 
to determine the Ni/Fe abundance ratios in multiple 
ejecta velocity components. We find consistent Ni/Fe 
abundance ratios that are ~5 times higher than the solar 
ratio, with an estimated uncertainty of 20%-—50%. Our 
reanalysis of Ni/Fe values from optical and near-IR data 
(which were previously interpreted as indicating much 
higher Ni/Fe abundance) using updated atomic data and 
allowing for local dust extinction shows that the revised 
values are roughly consistent with the JWST results. We 
compared the measured Ni/Fe abundance ratios with 
nucleosynthetic yields from Fe core-collapse and ECSN 
models available in the literature. We conclude that the 
modestly elevated Ni/Fe ratios in the Crab Nebula are 
consistent with either a low-mass Fe core-collapse 
explosion or an ECSN. The kick velocity of the Crab's 
pulsar more strongly favors an LMCCSN that can 
explode with equally low energy but a higher asymmetry 
than an ECSN. We noted that some models for low- 
energy Fe CCSNe at the lowest range of initial masses 
may provide a suitable mechanism for eruptive pre-SN 
mass loss that might lead to strong CSM interaction that 
could explain the high luminosity of SN 1054. 


Future studies combining this JWST data set with additional 
observations can include an exploration of the small-scale 
structure of the Rayleigh-Taylor filaments and mixing of 
ejecta. The photoionization models can be expanded to include 
the full set of lines detected with MRS to determine other 
elemental abundances that can be compared to explosion 
nucleosynthesis models. A more detailed study of the small- 
scale spatial and temporal variations in the synchrotron 
emission can be conducted with the NIRCam and MIRI 
images and compared with other multiwavelength data. 
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Additional spectroscopic observations sampling other locations 
within the ejecta filaments will be important for determining if 
the Ni/Fe abundance ratios vary across the filaments. Finally, a 
limiting factor in constraining the Crab Nebula’s origin is in the 
uncertainty of theoretical predictions, so refining and expand- 
ing theoretical explosion and nucleosynthesis models for 
progenitors at the low end of the mass range (8-10 M.) will 
be crucial for more definitively characterizing the progenitor 
and explosion that produced the Crab Nebula. 
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